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We compare mean-link and average plaquette tadpole renormalization schemes in the context of the quarkonium 
hyperfine splittings in lattice NRQCD. Simulations are done for the three quarkonium systems cc, be, and bb. The 
hyperfine splittings are computed both at leading and at next-to-leading order in the relativistic expansion. 
Results are obtained at a large number of lattice spacings. A number of features emerge, all of which favor 
tadpole renormalization using mean links. This includes much better scaling of the hyperfine splittings in the 
three quarkonium systems. We also find that relativistic corrections to the spin splittings are smaller with mean- 
link tadpoles, particularly for the cc and be systems. We also see signs of a breakdown in the NRQCD expansion 
when the bare quark mass falls below about one in lattice units (with the bare quark masses turning out to be 
much larger with mean-link tadpoles). 



Tadpole diagrams in lattice theories are in- 
duced by the nonlinear connection between the 
lattice link variables and the continuum gauge 
fields. This causes large radiative corrections to 
many quantities in lattice theories. Most of the 
effects of tadpoles can be removed by a mean field 
renormalization of the links H 



1t 



(1) 



where an operator dominated by short-distance 
fluctuations is used to determine uq. 

One of the earliest applications of tadpole im- 
provement was in the development of lattice non- 
relativistic quantum chromodynamics (NRQCD) 
j2| [7j. Precision simulations of the T system in 
NRQCD have provided important phenomenolog- 
ical results, including the strong coupling con- 
stant [Q and the 6-quark pole mass |||| . However 
the situation for charmonium is more problem- 
atic, due to large relativistic corrections |8). 

The quarkonium spectrum provides a powerful 
probe of tadpole renormalization. The quarko- 
nium hyperfine spin splittings in particular are 
very sensitive to the details of the NRQCD Hamil- 
tonian, with the relevant operators undergoing 
large tadpole renormalizations. For example, it 
has been shown || that scaling of the charmo- 
nium hyperfine splitting is significantly improved 



when the tadpole renormalization is determined 
using the mean-link uq l measured in Landau 
gauge 

u ,l = (ReTr , = 0, (2) 

compared to when the fourth root of the average 
plaquette uo,p is used: 

u q .p = <Wrf/ pl ^ 1/4 . (3) 

We make a comparison of the two tadpole 
renormalization schemes u .l and u p (imple- 
mented at tree- level) in the context of the quarko- 
nium hyperfine splittings in NRQCD. This is done 
for the three quarkonium systems cc, be, and 
bb. The hyperfine splittings are computed both 
at leading (0(Mqv 4 )) and at next-to-leading 
(0(Mqv 6 )) order in the relativistic expansion. 
(For further details see Ref. Q) 

All quantities are calculated after re-tuning of 
the lattice action parameters for each system. 
The resulting quark masses and lattice spacings 
for the three quarkonium systems for the NRQCD 
action at 0(v 6 ) are given in Tables [j] and |. To 
minimize systematic errors from quenching the 
6-quark mass is tuned separately to 1 Sq for the 
be and bb states with the c-quark mass tuned 
to the 1 So cc state. Likewise, the lattice spac- 
ings were determined separately from the spin- 
averaged IP — 15 mass difference, which we set to 
458 McV (experimental value for charmonium). 
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Table 1 

Bare quark masses for the three quarkonium sys- 
tems at 0(u 6 ), using Landau gauge mean- link 
tadpoles Uq,lI the stability parameter n for each 
mass is given in square brackets. 



(3l 


a C c (fm) 


aM c >] 


aM»[n] 


aM»[n] 


7.5 


0.155(4) 


1.10 [4] 


3.20[2] 


3.20[2] 


7.4 


0.179(2) 


1.20 [4] 


3.57[2] 


3.57[2] 


7.0 


0.280(4) 


1.97[2] 


6.10[2] 


5.35[2] 


6.85 


0.319(5) 


2.25[2] 


6.50[2] 


5.90[2] 


6.7 


0.361(6) 


2.50[2] 


7.20[2] 


6.35[2] 


6.6 


0.380(7) 


2.67[2] 


7.50[2] 


6.66[2] 


Table 2 



Bare quark masses for the quarkonium systems 
at 0(v 6 ) using average plaquette tadpoles i*o,p. 



pp 


a ce (fm) 


aM c u [n] 


aM»[n] 


aM»[n] 


7.3 


0.140(4) 


0.65[8] 


2.87[2] 


2.87[2] 


7.2 


0.169(2) 


0.83 [4] 


3.20[2] 


3.20[2] 


7.0 


0.210(2) 


1.10[4] 


4.10[2] 


3.95[2] 


6.8 


0.256(3) 


1.43 [3] 


4.98[2] 


4.53[2] 


6.6 


0.313(4) 


1.80 [3] 


5.83[2] 


5.23[2] 


6.4 


0.350(6) 


2.15[2] 


6.45[2] 


5.60[2] 


6.25 


0.390(6) 


2.41[2] 


6.85[2] 


5.99[2] 



The lattice NRQCD effective action for quarko- 
nium is organized according to an expansion in 
the mean squared velocity v 2 of the heavy quarks, 
with corrections included for lattice artifacts. 
The effective action, including spin-independent 
operators to 0(u 4 ), and spin-dependent interac- 
tions to 0(v 6 ), was derived in Ref. ||. Following 



Rcfs. 



we use the evolution equation 



G t+1 = 1- 



~2n 



uUi- 



aHc 
2n 



(l-aSH)Gt. (4) 



Relativistic corrections are organized in powers 
of the heavy quark velocity: 



5H = SH^+SH^ 6 \ 



(5) 



Only next-to-leading spin-dependent interactions 
are considered in §H( 6 > . 

Simulations were done with the derivative oper- 
ators and the clover fields corrected for their lead- 
ing discretization errors. Complete expressions 
for the operators can be found in Refs. P,Bl2j. 



Gauge-invariant source and sink smearing was 
used for Meson operators. 

The gauge-field configurations were made us- 
ing a tree-level (9(a 4 )-accurate tadpole-improved 
action f|l3ll 



Figure 1. Hyperfine splittings with iio,£ versus 
lattice spacing squared. 
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Figure 2. Hyperfine splittings with uo.p versus 
lattice spacing squared. 
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There are a number of clear features in the data 
as shown in Figures. |l| and||. To begin with, note 
that the results with uq.l show smaller scaling vi- 
olations than the results with uo,p ■ The small- 
est scaling violations are in uo.l at 0(v 6 ), which 
show little change over a large range of lattice 
spacing. The scaling analysis provides evidence 
that uq.l tadpole renormalization yields a more 
continuum- like action than does uq p. 
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The most striking feature is the drop in the be 
splitting at smaller lattice spacings, when uo,p is 
used at 0(v 6 ) . Most c-quark data with uo t p show 
large changes at small lattice spacings. The 
data exhibit much smoother behavior. 

We interpret these features as possible indica- 
tors of a breakdown in the NRQCD effective ac- 
tion at smaller lattice spacings, when the bare 
quark mass in lattice units ciMq falls below one. 
The bare c-quark mass is larger when mo,l is used. 

Another key feature is that the relativistic cor- 
rections to the hyperfine splittings are smaller 
when the action is renormalized using uq.l- For 
example, we find that the charmonium hyperfine 
splitting is reduced by about 30-40% in going 
from 0(t> 4 ) to 0(v 6 ) when using «o,l, compared 
to a reduction of about 40-60% when using Uo,p- 
(Relativistic corrections have been analyzed in 
the T system |T^,|l5[, and in heavy-light mesons 
|l6|-^8[). The u ,p relativistic corrections depend 
strongly on the lattice spacing and may be re- 
lated to the pathologies discussed above. This 
casts new light on the results obtained in Ref . || , 
where relativistic corrections to spin splittings in 
NRQCD were first calculated. The velocity ex- 
pansion for charmonium may not be as unreliable 
as was suggested. 

We note finally that it is reasonable to attempt 
to extrapolate the 0(v e ) hyperfine splittings for 
cc and be to zero lattice spacing, from the data 
on coarse lattices, where there is good scaling be- 
havior. However, the extrapolations in the Uq,l 
and uo.p data are clearly very different. This sug- 
gests that some relevant operator coefficients a in 
the NRQCD action receive significant 0(a s ) cor- 
rections. This underlines the need to go beyond 
tree-level tadpole improvement in order to clarify 
the differences between renormalization schemes. 

We have presented evidence that favors tadpole 
renormalization using the mean-link in Landau 
gauge over the fourth root of the average plaque- 
tte. This includes a demonstration of better scal- 
ing behavior of the hyperfine splittings in three 
quarkonium systems when uqx is used, and a 
smaller size for spin-dependent relativistic correc- 
tions. These results help to elucidate the struc- 
ture of the NRQCD effective action. 
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